Sticky Surfaces: Adhesion Based Microbial Aggregation Dynamics 
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We consider spherical particles coated with binding ligands and study their attachment in quiescent 
flow. This fluid- immersed adhesion is widespread in many biological and artificial settings. Our 
theory highlights how the physics of the binding kinetics of these ligands (expressed through the 
collision factor function) as well as the attractive / repulsive surface potential in an ionic medium 
effects the eventual size of these particle aggregates (floes) . For experimentally measurable parame- 
ters in microbial population studies, our results suggests that large floe aggregates are possible with 
more stretchable ligands due to efficient inter- floe collisions (or a large, non-zero collision factor). 
Strong electrolytic composition of the surrounding fluid favors large floe formation as well. 

PACS numbers: 



Introduction - The formation of aggregates, induced by 
the adhesion of two spherical particles or nearby surfaces 
is important in many scientific and industrial processes. 
In particular, interfacial attachment leading to larger floe 
aggregates via the latching of binders on surfaces in close 
proximity is widespread. Examples include binding of 
bacterial clusters to medical implants and host cell sur- 
faces [1], cancer cell metastasis [2], and the coalescence 
of medical gels with nano-particles for targeted drug de- 
livery [3]. Moreover, coagulation and flocculation (the 
chemical and the physical aspects of adhesion) are also 
important in pulp and paper-making industries as well as 
waste water treatment plants [4j|. Although the micro- 
scopic structure and geometry of the adhering surfaces 
of floes in these phenomena vary greatly, they share a 
common underlying physics. Past investigations in theo- 
retical modeling of fluid-borne surface adhesion have ad- 
dressed ligand-receptor binding kinetics jjj Q , surface de- 
formation 0, [H and flow past the surrounding surfaces 
d, [To| . Our aim in this letter is to explore how this ad- 
hesion {collision in colloid science literature) mechanism 
for spherical floes is governed by various geometric and 
fluid parameters as well as how the surface forces and 
binding kinetics of the ligands impact the eventual sizes 
of these floes. 

Model - Our future studies are geared towards tracking 
bacterial aggregates in suspension, as a continuum mass 
of EPS (extracellular polymeric substance) network 
Following the general outline given in [12|, we define 



b(t, x, s) Ax =number of bacterial biomasses having volu- 
mes between x & x + Ax in time t at a sp- 
atial point s. 

In volumes between x\ and X2, the total number of 
floes Bq is given by 



for [x±, X2] C [x, x], where x and x are the minimum and 
maximum aggregate volume sizes, respectively. A finite 
nutrient supply and the duration of the experiment allow 
us to assume that x is finite. Further, the extracellu- 
lar structures must be actively maintained and thus the 
minimal size x is the volume of one bacteria. The conser- 
vation of bacterial biomass through aggregation, or the 
governing equation for b is [12| 



b t + v • V s b = A(x, b) 



(2) 



px 2 

Bo(t, xi, X2, s) = / b(t,x,s)dx 



(i) 



where v is the relative velocity of one floe with respect 
to the other. A = Ai n - A out , where A- in is the rate with 
which floes of size in [x, x + Ax] are created and A out is 
the rate a floe of size in [x, x + Ax] joins with another 
floe, to form a volume greater than x + Ax. These rates 
are given by 

2 r-X — X 

A in (x, b) = - J K A (y, x - y)b(t, y, s)b(t, x-y, s)dy, 

x G [2x,x] (3a) 
AmtO, V) = b(t, y, s) / K A (x, y)b(t, y, s)dy, 

J x_ 

x G [x, x — x] (3b) 

Ka is the aggregation kernel, describing the rate with 
which floes of volume x and y combine to form a floe of 
volume x + y. The next two sections will focus on model- 
ing this kernel based on the hydrodynamics and surface 
potentials of two coalescing spherical floc-surfaces. 
Adhesion mechanics - We present a simple model of inter- 
facial attachment between two spheres (of radius Ri and 
R2, Fig JT]) immersed in a fluid medium [H[. The surface 
of the spheres bind onto each other due to the presence 
of adherent elastic binders (polymer strands with sticky 
heads) on the surfaces, as well as the attractive surface 
potential. The adherent binders are idealized as linear 
Hookean springs with stiffness ^0 and mean rest length 
lo (FigHJ). The effect of the shear flow on the mean rest 



2 



length of the binders, as well as the shearing effect of the 
mean flow on the binders is neglected [HI . These simpli- 
fications allow us to focus on the dynamics normal to the 
surface, but could be relaxed for a more general caseQ. 



due to the surface potential. The direction of this force 
is along the direction vector from the spherical floe of 
radius Ri to the floe of radius R2. The total force arising 
from all such bonds is given by 
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FIG. 1: Two polymer coated coalescing spherical floes. 

For a given spatial point s = (si, 52,53), define 
ATot#(s, t)dA as the number of bonds in the transverse 
direction that are attached between the surfaces dA at 
time D(s) be the minimum distance between the two 
spheres, A^ot be the total number of binding ligands and 
g be the density of bound ligands on the adhesion surface. 
In floe literature, the function g is synonymous to the 
collision factor as well. Hence the total number of bonds 
formed in the transverse direction is J A ATot#(s, t)cL4, 
with A c being the area of adhesion. The bond attache- 
ment/ detachment rates, influenced by the surface poten- 
tial of the two charged surfaces, are 



K on (s) = iC n exp [ 
Koff(s) = iCff ex P [ 



-K s (D(s)-l ) 2 ^W(D(s)) 

2k B T J 
(acq - K 8 )(D(s) - h) 2 + W(D(s)) - 
2k B T 



(4) 



where k B is the Boltzmann constant, T is the tempera- 
ture, k s is the spring constant of the transition state used 
to distinguish catch (k < k 3 ) from slip (k > k 3 ) bonds 
[Ej], W(D) is the total surface potential, if* n eq ,i^* ff eq 
are the equilibrium binding affinities. For notational sim- 
plicity, we denote D(s) = D. In the limit of small binding 
affinity and abundant binding receptors on the surface of 
adhesion (^Tot-^on, ^/K oii} eq <C 1), the bond ligand density 
evolves as the following differential equation (EJ: 



dt 



^Tot^on 



K sg, 



(5) 



where = ^| + v • V s g. The instantaneous force that 
these two colliding charged surfaces exert on each other 
(and acting normal to the surface) becomes: 



f(M) = «o(2?-fc) + V D -W%D) 



(6) 



where the first term represents the stretching force due 
to Hooke's law and the second term represents the forces 



F(s, t) = A Tot [ g(s, t)f (s, t)dA(s, t) (7) 

In the limit of small binding affinity and normal adhesion, 
the adhesion area, A c , is given by A c = tt R 2 [13[, where 
the adhesion radius, R C1 is 

R c (s,t)= (^) 1/ \(i?- 1 / 2 (s,t) + i? 2 - 1/2 (s,0) (8) 
Finally, in a Stokes flow, the aggregation rate, K^, is 



K A = 7aA c f/ c 



(9) 



( is the drag coefficient, 7^ is the aggregation contact effi- 
ciency parameter. Eqns. ([3a] \3b\ \2\ \5\ l6| \7\ l9j) along with 
initial conditions, 6(0, x) = bo(x), is the entire system. 
Long range interactions - Experiments are being per- 
formed to determine the physiochemical characteristics 
of two coalescing, polymer-coated, charged bacterial floes 
immersed in ionic solvents [14] . We describe these inter- 
actions through the DLVO approach, i.e., the Coulom- 
bic and Van der Waals interaction. To simplify modeling 
framework, other interactions including hydration effects, 
hydrophobic attraction, short r ang e steric repulsion and 
polymer bridging are neglected [15j . 

For two charged spheres, with radii Ri, R2, the repul- 
sive Coulombic forces in the gap, D, is given by 



W C (D) = 27re e^k B T 



/ 2R 1 R 2 - 



(10) 



where k is the Debye length, e, eo are the dielectric con- 
stant of vaccum and the medium, respectively, ^0 is the 
zeta potential or the electric potential of the diffuse cloud 
of charged counter ions. The attractive Van der Waal 
forces for spherical floes in the regime of close contact 
#2), is 



W VW (D) = 



Ak B T R X R 2 
6D ft! R 2 



(11) 



where A is the Hamaker constant, measuring the van der 
Waal 'two-body' pair-interaction for macroscopic objects. 
The net surface potential is W(D) = W C (D) + W VW (D). 
This potential is pair- wise attractive over very short and 
long distances, and pair- wise repulsive over intermediate 
distances (Fig|2]). 

Results - We consider the case of homogeneous floe aggre- 
gation (the variables are independent of the spatial loca- 
tion, s), under quiescent flow conditions. If the adhesion- 
detachment rate of the floes is sufficiently rapid so that 
they remain in equilibrium, the fluid surrounding the 
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floes is stagnant and the effects of peeling velocity is neg- 
ligible [5], then the bond-ligand density (i.e. the collision 
factor), g, (Eqn. SJ[5]) evolves according to 
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iTot 



K 



^° 2fc R T 



K* 



(12) 



l off "off 

and the total adhesion force between two floes reduces to 
F(t) = A Tot gt[7rR 2 c (t)} (13) 
The parameters used in simulations are listed in Table U 



Parameter 


Value 


Units 


Source 




(0.01-10) xlO" 5 


N m -1 


M 


fa 


io- 8 


m 


[13] 


c 


(0.01-2.5) 


N m" 1 s 




A Tot 


10 y 


m~' 2 


[ia 


K on/K* ff 


io- 12 




[13] 




2.7 xlO -1 ^ 


fL- 2 


[12] 


eo 


8.854 xlO" 12 


Fm- 1 




e 


78.5 






^0 


25 


mV 


[15] 


ft 


1.5 


nm" 1 


[15] 


A 


2.0 




[15] 



TABLE I: Parameters common to all simulations. 

Surface potential and collision factor - The pair-wise 
surface potential, W(D), is valid over short distances 
(D <C R). The electrolyte is assumed to be a 1-1 elec- 
trolyte at a concentration of 50 mM, zeta potential of 
25mV and Debye length for screened charges n~ l = 1.5 
nm -1 , a typical of values for particles in water (l6| . 

With increasing binder stiffness, the potential energy 
barrier as well as the depth of the secondary minima 
increases (Fig. [2^l) . Hence aggregation becomes less fa- 
vorable with increasing binder stiffness, a fact seen in 
FigHK, as well. However, at critical separation distances 
(5 < D < 12 nm), the depth of the secondary minima 
remains lower than the average thermal energy of the 
particles (3fc#T/2), implying that at this range the ag- 
gregating floes are colloidally stable. The primary min- 
ima shown here is unphysical, since at very short separa- 
tion distances the non-DLVO repulsive steric interaction 
is dominant that prevents the particles from coming into 
true contact. The regions of attraction/repulsion of this 
potential is inferred from surface force, f (Fig. Wp>)- At 
intermediate distances (lnm < D < fa), the repulsive 
Coulombic forces are dominant while at long distances 
(D > fa), the adhesive forces are dictated by the attrac- 
tive spring force of the stretched binders. We choose 
D* = llnm, where the adhesive forces are attractive. 

The collision factor, g(D), (where the separation dis- 
tance, D = Sl ft S2 + fa, with the center of the spatial frame, 
fixed at the point of minimum separation on sphere 2, R 
being the radius of identical colliding floes, see FigH]) is 
symmetric about the mean rest length of the binders, fa, 
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FIG. 2: (a) Total pairwise potential, W(D), versus the floe 
separation distance, D, for two identical spherical microbial 
floes with varying binder stiffness but fixed size, R=l/mi, (b) 
Surface potential force, f, for different binder stiffness, show- 
ing regions of attraction (f > 0) and repulsion (f < 0) 



in the absence of fluid flow (Fig|3]). Since the floe separa- 
tion distance with a significant non-zero contact is consid- 
erably reduced for stiff binders, the adhesion-detachment 
mechanism is more efficient for tensile springs. 
Numerics - To solve our model, we employed the dis- 
cretization scheme developed by Banks [17j and adopted 
by Doumic (T^j. The convergence of the scheme was 
tested using the test functions in [l2j|. A linear rela- 
tionship between the L°°-error and the mesh-size, Sx 
was found using this first order approximation scheme. 
The initial number density is chosen as bo(x) = 7.47 x 
10 _4 e _0 00676:c , where the coefficients are fit to the exper- 
imental data from the Younger Lab [12j . All solutions are 
marched in time until T=100 minutes. We chose 1 fem- 
toliters (fL) as a lower bound x in our simulations. Our 
aggregation model allows the upper bound, x, of the do- 
main to go unrestrained (i.e., x — >• oo), but the results 
are presented inside the window 1 < x < 1000 fL. 

In this letter we limit our parameter analysis to binder 
stiffness and the electrolytic composition of the fluid 
medium. Fig. [4^ suggests that stiff binders lead to fewer 
large aggregates (i.e., b(x, T, = 10 -2 ) < b(x, T, 
= 10" 3 ) < b(x, T, ft = 1(T 5 ), for x > 600). Since 
aggregation is influenced by the collision factor, a higher 
value of g(D*) suggests that two floes close to each other 
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FIG. 3: Collision factor g vs spatial coordinates (si,S2) for 
identical colliding floes of radius K—lum with (a) ko = 10 -2 
Nm" 1 and (b) = 10 -5 Nm" 1 . Stiff ligands have a nar- 
row range of separation distance D(si,S2) with a non-zero 
collision impact. 
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FIG. 4: Floe number density distribution versus floe- volume 
at (a) different binder stiffness and screening length, k = 1.5, 
and (c) different screening lengths and no = 10 ~ 3 Nm -1 . The 
dash-dot curve in these figures is the initial conditions (bo(x)). 

are more likely to coalesce (Fig. [3]). This factor is greatly 



compromised at a floe separation distance, D* = llnm, 
for stiff binders, leading to fewer large aggregates. Simi- 
larly, the number of large aggregates drops for a weaker 
electrolytic fluid composition (or longer Debye length, ft, 
Fig.Ht>), suggesting that a weaker ionic composition does 
not favor aggregation either. 

Conclusion - We have presented a minimal aggrega- 
tion model of coalescing spherical floes under the influ- 
ence of surface potential. Attempts are made to explain 
the physical implications of the 'collision factor' a term 
widely used but scarcely understood in engineering lit- 
erature |4[. Predictions about the floe aggregate size, at 
various fluid and surface potential parameters, have been 
made using numerical simulations. These predictions are 
being currently validated with a series of experiments, 
e.g., those detailing changes in the floe sizes at various 
ionic composition of the static fluid, or with different 
polymer coated floe surfaces. 
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